Strong surface scattering in ultrahigh mobility Bi 2 Se 3 topological insulator crystals 
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While evidence of a topologically nontrivial surface state has been identified in surface-sensitive 
measurements of Bi2Se3, a significant experimental concern is that no signatures have been observed 
in bulk transport. In a search for such states, nominally undoped single crystals of Bi2Se3 with 
carrier densities approaching 10 16 cm -3 and very high mobilities exceeding 2 m 2 V -1 s _1 have 
been studied. A comprehensive analysis of Shubnikov de Haas oscillations, Hall effect, and optical 
reflectivity indicates that the measured electrical transport can be attributed solely to bulk states, 
even at 50 mK at low Landau level filling factor, and in the quantum limit. The absence of a 
significant surface contribution to bulk conduction demonstrates that even in very clean samples, 
the surface mobility is lower than that of the bulk, despite its topological protection. 

PACS numbers: 71.18.+y, 72.20.My, 78.30.-j, 78.20.Ci 



Topological insulators are bulk insulators that feature 
chiral Dirac cones on their surface. This surface state, 
which is protected by time-reversal symmetry, is of fun- 
damental interest and may have potential for spintron- 
ics and quantum computation applications. A small 
group of semimetallic stoichiometric chalcogenides have 
been theoretically shown to have the properties of three- 
dimensional topological insulators: Bi2Se3, Sb2Te3, and 
Bi2Te3 [HQ. The chiral surface state has been investi- 
gated by surface probes, notably angle-resolved photoe- 
mission spectroscopy (ARPES) 2-4] and scanning tun- 
neling microscopy (STM) la, [6J and similar studies exist 
of the alloy Bii-^Sb^ [7|, |8|. The observation of chiral 
Dirac cones and forbidden backscattering in these exper- 
iments has generated a great deal of excitement. 

In principle, electrical conduction in topological insula- 
tors should occur only at the surface, but in practice these 
stoichiometric materials have been known for decades to 
be low carrier density metals. Bi2Se3 is expected to have 
a 300 meV direct bandgap at zone center [l|, Q , yet is al- 
most always n-type, long believed the result of charged Se 
vacancies. It remains an open question whether a bulk in- 
sulating state can be achieved in stoichiometric undoped 
Bi 2 Se3. This problem has been recently sidestepped by 
counter-doping with Ca @, H, [l(| , although this proce- 
dure introduces further defects. Over a range of carrier 
densities, ARPES data show one conduction band at zone 
center 0, IH > and bulk Shubnikov de Haas measurements 
indicate that the Fermi surface is ellipsoidal [llj . 

There have been several recent studies of bulk trans- 
port in Bi 2 Se3. A thorough study of angle dependent 
Shubnikov de Haas (SdH) oscillations on a sample with 
a moderate carrier density has mapped the anisotropy 
of the Fermi surface A comparative investigation 

demonstrated that while ARPES may identify a Fermi 
level in the band gap, SdH oscillations in similar samples 
always show metallic behavior, suggesting that the Fermi 
level may vary between bulk and surface in Bi2Se3 [l3j . 



This raises questions about the correspondence between 
surface and bulk studies on the same material. While it 
is possible to tune the Fermi level of Bi2Se3 into the gap 
via Ca doping, the resulting samples show unusual mag- 
netoresistance fluctuations, but no SdH oscillations from 
the surface state [ltj. In contrast, field-dependent oscilla- 
tions ascribed to the surface state have been observed in 
tunneling spectra measured on Bi2Se3 thin films 3, HI 
and Aharonov-Bohm interferometry of the surface state 
is reported in Bi2Se3 nanoribbons [l6|. Because SdH os- 
cillations attributable to surface states were identified in 
intrinsically disordered Bii_ x Sba; alloys [l7| . it is intrigu- 
ing that analogous bulk transport phenomena have not 
been seen in stoichiometric Bi2Se3. 

In order to look for signatures of the novel surface 
state in bulk nominally undoped Bi2Se3, single crystals 
were synthesized over a wide range of carrier concentra- 
tions. Their longitudinal and transverse electrical resis- 
tivity were systematically characterized as a function of 
field and temperature, and infrared reflection and trans- 
mission were studied. This investigation confirms that 
all measured transport properties can be ascribed to bulk 
electron-like carriers, despite the samples having the low- 
est reported bulk carrier densities and highest mobili- 
ties. Even at low Landau level filling factor and in the 
quantum limit, no transport signature of surface states 
is apparent, placing strong constraints on the mobility of 
surface carriers in these high quality samples. 

Single crystals of Bi2Se3 were prepared by melting 
high purity bismuth (6N) and selenium (5N) in sealed 
quartz ampoules. Multiple batches were prepared with 
varying bismuth/selenium ratios and heating conditions, 
which were responsible for variations in the carrier con- 
centrations of the samples. Similar trends were reported 
decades ago [l8[ and more recently The resul- 

tant crystals cleaved easily perpendicular to the C3 trig- 
onal axis. Four-probe measurements of longitudinal and 
transverse electrical transport between 1.8 K and 300 K 
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were conducted using a Quantum Design 14 T PPMS. 
The current flowed in the plane perpendicular to C3 in 
both longitudinal and transverse geometries, and the 
magnetic field H was always applied perpendicular to 
applied current. Electrical transport measurements at 
lower temperatures were performed in a dilution refriger- 
ator with a rotating sample stage and 15/17 T magnet. 
Optical measurements were performed using a Bomem 
DA3 FTIR spectrometer. Carrier density n was deter- 
mined by measurements of Hall effect, SdH oscillations, 
and fits to optical reflectivity data, which always yielded 
identical results within experimental uncertainty. 

The temperature dependence of the longitudinal elec- 
trical resistivity p(T) is summarized in Figure QJi. The 
values of n in these samples were estimated via analy- 
sis of the SdH oscillations in p{H) and approximating 
the Fermi surface as spherical, which is sufficient for pur- 
poses of comparison at low n [ll[ . In sample vi, no SdH 
oscillations are observed up to 14 T, and the estimated 
n = 5 x 10 16 cm~ 3 is based on Hall measurements on sam- 
ples from the same batch (Fig. [TJd). For n > 10 18 cm~ 3 , 
the p{T) data reflect good metallic conductivity and gen- 
erally, as n decreases, the samples become more resistive. 
For n < 10 17 cm -3 , a negative slope develops above 
250 K, which reflects a crossover to activated behavior 
at T > 300 K 19]. For n < 10 18 cm" 3 , the p(T) data 
also develop a shallow local minimum at 30 K, which is 
associated with a relatively small upturn that saturates 
as T -> 0. 

In Figure [TJ), the T dependence of the Hall carrier den- 
sities iih — (i?ne) _1 is shown for a range of samples. 
Values of the Hall coefficient Rn were determined via 
linear fits to symmetrized transverse magnetoresistance 
data, which are linear to 5 T over the entire T range, 
except for at lowest n where linearity extends to 2 T at 
low T. Generally, nn(T) has a gentle T-dependence, con- 
sistent with a gradual crossover from extrinsic to intrin- 
sic conduction (E g > 100 meV) 19]. However, metallic 
p(T) (Fig.QJ,) is always observed, indicating that phonon 
scattering dominates any increase in conductivity due 
to increasing carrier number for most of the T range. 
A detailed comparison of p(T) (sample Hi) and nn(T) 
between two samples from the same batch is shown in 
Fig. [2K These data highlight two regimes: below 30 K, 
where njj is constant and p(T) exhibits an upturn, and 
between 150 and 250 K, where there is a change in cur- 
vature in nn(T 1 ) and p(T). The magnitude of these lat- 
ter features is greatest in samples with n < 10 17 cm" 3 . 
The two T ranges are accentuated in a plot of the Hall 
mobility p(T) = R H (T)/p(T) (Fig. [5b), highlighting the 
substantial low-T mobility p > 2 m 2 V -1 s _1 in these 
samples, which is the highest value reported for Bi2Se3 
and corresponds to a mean free path longer than 300 ran. 
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A comparison to p(T) of batches with n ~ 10 
and n ~ 10 19 cm -3 indicates that p does not simply in- 
crease with decreasing n, although at low n it still exceeds 
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FIG. 1. (Color online) a) Comparison of the electrical re- 
sistivity p(T) between samples of Bi2Se3. Carrier densities 
n (cm -3 ) are estimated from SdH oscillations: i 1 x 10 19 , 
ii 5.3xl0 18 , Hi 4.9xl0 17 , iv 3.7xl0 17 , v 3.3xl0 17 , vi ~ 10 16 . 
In samples with n < 10 18 cm -3 , a shallow minimum develops 
at 30 K. b) Comparison of Hall carrier density riH between 
various samples. Low temperature values span 3 orders of 
magnitude. 



1 m 2 V- 1 S" 1 . 

In Fig. p(H) of sample Hi is shown. In addition 
to the dramatic anisotropy, oscillations are readily ap- 
parent. The non-oscillatory part is sufficiently well de- 
scribed by quadratic polynomial or weak power law func- 
tions. In contrast, sample vi exhibits no oscillations. 
The T-dependence of SdH oscillations is shown in Fig. [3] 
in panels a) and b) for two different field orientations 
(c.f. Fig. [2J;). The SdH oscillations have a frequency 
JsdH = 20 T when H \\ C3 and /sdH = 25 T when 
H _L C3, which correspond to an ellipsoidal Fermi sur- 
face with n = 6.3 x 10 17 cm -3 . A Lifshitz-Kosevich 
fit to the T-dependence of oscillation amplitude yields 
masses mil = 0.14(l)m c and m± = 0.16(l)m o , where 
m e is the bare electron mass, consistent with other 
12, [lij ]. In order to check whether lowering 
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the temperature would uncover a second frequency of 
SdH oscillations arising from the surface state, measure- 
ments were performed at 50 mK on two samples with 
riH = 4 — 5 x 10 17 cm~ 3 from the same batch as sample 
iv. The field angle dependence of both longitudinal and 
transverse resistance are shown in Fig. [3] panels c) and 
d). Figure [3J; shows oscillations in the transverse magne- 
toresistance R xy (H), whose frequencies range from 20 T 
with H || c 3 to 27 T with H 1 c 3 (n = 7 x 10 17 cm- 3 ). 
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FIG. 2. (Color online) a) A comparison of electrical resistivity 
p(T) (Hi from Fig. [1]) and Hall carrier density nn(T) in two 
samples from the same batch. The minimum in p(T) coincides 
with the saturation of ?ih(T) below 30 K. Despite increasing 
riH, p increases with T up to 300 K. b) The calculated mobility 
exceeds 2 m 2 V" 1 s" 1 at low T (labeled 10 17 ). The mobilities 
of samples with lower and higher n are compared, c) Field 
orientation dependence of p at 1.8 K. Compare to p(H) of 
sample vi (rescaled by factor of 10), which exhibits no SdH 
oscillations. 



FIG. 3. (Color online) Temperature dependence of SdH os- 
cillations with a) H || C3 trigonal axis and b) H _L C3. Os- 
cillation frequencies /sdH = 20 T and 25 T, respectively, and 
the mass msdH ~ 0.15m c . SdH oscillations from two differ- 
ent samples at 50 mK in the c) transverse and d) longitudinal 
MR. The largest amplitude and smallest frequency in c) is 
observed with H || C3 and the frequency ranges from 20 T 
to 27 T. In d) the largest amplitude is observed with H _L C3 
and /sdH = 15 T. Only one SdH frequency is observed in each 
measurement. 



As expected, the magnitude of the effective Hall signal 
diminishes as H rotates into the plane and the evolution 
of the oscillation frequency reflects the expected Fermi 
surface anisotropy. The longitudinal SdH oscillations are 
shown in Fig. [3J1. The most prominent oscillations, with 
/ sdH = 15 T (n = 3 x 10 17 cm" 3 ) occur with H _L c 3 . 
In this field orientation, the electrons are almost all in 
the first Landau level. Weak oscillations with H || C3 
have /sdH = 11 T, yielding an anisotropy close to that 
of the other two samples. The lack of bulk SdH oscilla- 
tions in sample vi (Fig. [2J;) having even lower n is due 
to its low estimated /sdH = 4 T, meaning that most of 
the accessible field range is in the bulk quantum limit. 

Samples were also characterized by infrared reflection 
and transmission, which confirm the estimates of n based 
on Hall and SdH data. Figure shows an experimental 
reflectivity spectrum of a sample from the same batch as 
sample iv. Data were taken at 6 K and are presented 
with a fit to a sum of Lorentz oscillators to model the 
complex dielectric function. The fit shown in Fig. |4^, in- 
cludes only bulk carriers as a Drude-type term together 
with two low frequency phonons. A strong low frequency 
phonon is centered at 67 cm -1 with a plasma frequency 
of 640 cm -1 and relaxation rate of 5 cm -1 and a sec- 



ond weaker phonon is centered second at 124 cm -1 with 
a plasma frequency of 76 cm -1 and a width of 2 cm -1 . 
The calculated optical conductivity is shown in Fig. 2b 
together with two dc values (sample iv) for rough com- 
parison. This yields a bare plasma frequency of free carri- 
ers to be 382 cm' 1 , corresponding to n = 2.4 x 10 17 cm~ 3 
for 0.15m c . The effective scattering rate 7 = 8 cm -1 
(r = 0.8 ps) is unusually low for bulk carriers, consistent 
with a high mobility sample. Figure 0fc shows the real 
part of the dielectric function. The fit parameters are in 
agreement with those from the transmission spectra of 
very thin crystals. An important difference, however, is 
that in transmission, a very broad Drude surface contri- 
bution (7 = 100 cm -1 ) is detected. 

For small n, the competition between negative E\ of the 
free carriers and the positive £1 from the phonons gives 
rise to an additional low frequency plasma edge {e\ = 0) 
that leads to a transmission window near 30 cm -1 . The 
strong bismuth-dominated low frequency phonon pro- 
duces a large static dielectric constant £i(0) s» 100 for 
Bi 2 Se3. This implies a large reduction in the Coulomb 
potential and thus the scattering rate from ionized impu- 
rities, which can account for the small scattering rate of 
the bulk carriers and the good metallic behavior of this 
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FIG. 4. (Color online) Spectra of a) reflectivity, b) optical 
conductivity, and c) dielectric constant. Curves b) and c) 
were calculated from the fit to the reflectivity data in a). For 
comparison, the blue circle and red square in b) are dc values 
from sample iv for 6 K and room temperature, respectively. 

doped semiconducting material, and protecting against 
localization at low n. The Coulomb mediated electron- 
electron interaction is also reduced [20| and, along with 
the presence of low frequency optical phonons, implies 
a dominance of the electron-phonon interaction both for 
the carrier mobility and for the electron-electron interac- 
tion. 

ARPES measurements show that surface states in 
Bi2Se3 always have a larger Fermi surface than the bulk 
and that the areal difference increases as the chemical po- 
tential decreases Quantum oscillations arising from 
the surface should have an easily identifiable higher fre- 
quency of about 300 T (fcp = 0.1 A -1 ) when the chemi- 
cal potential is near the bottom of the conduction band. 
However, no higher frequency oscillation is ever observed 
in our measurements, even in the quantum limit, where 
no bulk SdH oscillations exist to obscure surface signal. 
As the existence of the surface states is not in doubt, the 
absence of their transport signature can only arise from 
very high scattering rates. A conservative assumption 
that the chemical potential is the same at the surface 
and bulk, consistent with ARPES surface aging data [ij , 
yields n sur f » 10 13 cm" 2 . In our samples, the bulk areal 
carrier density is about 100 times that of the surface, 
while the peak-to-noise ratio in the Fourier transform of 
the SdH oscillations is greater than 500, so the surface 
mobility has to be at least 5 times worse than the bulk 
mobility for SdH oscillations to be below our detection 
threshold. We calculate scattering times Tbuik = 2 ps and 
Tsurf < 0.4 ps. The latter value is generous compared to 
an upper limit of r su rf < 0.05 ps determined from our 
transmission measurements. These estimates indeed in- 
dicate substantial surface scattering, yet are entirely con- 



sistent with ARPES data, which show significant surface 
band broadening, with t = 0.02 ps Q or 0.04 ps [2l[ even 
for samples cleaved in high vacuum. 

It is clear that topological protection does not guaran- 
tee high surface mobility, so increasing the surface con- 
duction will require considerable effort. Reducing sample 
thickness is one possibility, but bulk transport proper- 
ties appear to be insensitive to sample thickness down 
to 5 (im [l3j]. Even in 50 nm thick nanoribbons, the ob- 
served / sdH = 70 T and 110 T are well below 300 T 
and are not likely due to surface states. Relative sur- 
face conduction may alternatively be improved by fur- 
ther reducing bulk n. Attempts to do so via chemical 
doping yield no surface SdH oscillations [l(|, possibly 
due to even further decreased mobility. On the other 
hand, doping may improve the surface signal if it pref- 
erentially decreases the bulk mobility without reducing 
surface mobility, which may be possible considering the 
already high surface scattering rate. A cleaner approach 
is to instead preferentially improve the surface mobility 
via controlled surface preparation. Such a task will be 
challenging, as it is difficult to improve upon the vac- 
uum cleaving employed in ARPES studies, which yields 
broad surface state linewidths. Nonetheless, high surface 
mobility is a prerequisite for precision experiments like 
optical Kerr rotation and future applications potential of 
topological insulators. 

We note a recent report of the observation of surface 
SdH oscillations in pulsed fields [22j |. The cleaved sur- 
faces of the Sb-doped Bi2Se3 were minimally exposed, yet 
the requirement of high magnetic fields points to strong 
surface scattering even under controlled conditions. 
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